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The ionic conductivity of mixed alkali glasses exhibits a deep minimum as a function of the relative 
concentrations of the two alkali ions. To study this behaviour we consider a simple one-dimensional 
model for asymmetric diffusion of two kinds of particles. Different particles are assumed to repulse 
each other. We consider two versions of the model: with or without overtaking of particles. For 
the case of perfect repulsion we find exact expressions for the stationary current. The model with 
weaker repulsion is studied by means of numerical simulations. The stationary current as a function 
of the ratio of particle concentrations is found to exhibit a minimum, related to correlations existing 
in this system. 
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I. INTRODUCTION 



Interesting physical effects have been observed in studying the behaviour of glassy ionic 
conductors. The conductivity of mixed alkali systems exhibits a strong dependence on the 
ratio of the alkali concentrations with a deep minimum at a ratio close to 1 ffl]. In some 
experiments the minimum conductivity is a few orders of magnitude smaller than the con- 
ductivity of the pure system. An interesting model explaining the anomalous conductivity 
has recently been introduced in (see also ||, 0]). The model is based on the assumption 
that cations in glasses create and maintain their own local environment which, due to mem- 
ory effects, produce a strong dependence of the conductivity on the relative concentrations 
of the two ions. Numerical studies of this model have been presented. 

Problems of this kind can be understood within the framework of lattice gas models. 
Recently an extensive effort has been invested in studying particles hopping in a preferred 
direction with stochastic dynamics and hard core interactions. These are simple exam- 
ples of non-equilibrium macroscopic systems || , || which exhibit very interesting collective 
phenomena such as phase transitions. These models have been used to study hopping con- 
ductivity or diffusion in narrow pores ||, |7|. Furthermore, they are related to growth 
processes ||- |]IT| . Some exact results for particular models are known ||- pl| . 

The aim of this paper is to discuss a lattice gas model for asymmetric particle hopping 
which is related to the behaviour of mixed ionic conductors. The model describes the 
dynamics of two kinds of particles moving in the same direction under the influence of a 
driving field. Here, we do not consider site memory effects but rather assume that there 
exists a repulsive interaction between the two types of ions. Our model, studied in d — 1 
dimensions, is simpler than the one introduced in 0. 

We consider two versions of the model. In the first version the two kinds of particles 
jump independently to the right along a one- dimensional lattice without overtaking each 
other. In the second version different kinds of particles are allowed to overtake. This may 
incorporate the effect of the higher space dimension into the one-dimensional model. 
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II. THE MODEL 



To model the dynamics of a system with mixed charge carriers we consider a one- 
dimensional lattice of length L. Each site can take one of three states: it may be occupied 
by a particle of type a or by a particle of type b, or it may be vacant (occupied by a hole e). 
We assume that both kinds of particles contribute to the conductivity of the system. The 
a- and b- particles may be identified with different alkali ions. The particles are assumed to 
undergo an asymmetric exclusion dynamics. The asymmetry is caused by a strong electric 
field. Throughout this paper we assume periodic boundary conditions. 

The total number of a- (b-) particles is assumed to be equal to p a L ( p b L ). We define r 
as the ratio of concentration of a-particles to the sum of both particle concentration. 

r = p a /p with p = Pa + Pb . (2.1) 

To be specific we consider stochastic asymmetric exclusion dynamics with interaction 
between the a- and 6-particles. Two versions of the model are studied. In the first version 
no overtaking of a- and 6-particles is allowed. A particle can move to its right if the site on 
its right is vacant. The probability of making the step during the time interval dt depends 
on whether an afe-bond is broken (adt), created ((3dt), neither broken nor created (ldt) or 
one bond is broken and another one created (afldt). The interaction between the two types 
of particles is thus introduced by the two parameters a and (3. The possible steps defining 
the dynamics are given by the following processes: 



xaex ==>- xeax with rate 1 (2.2) 

baex =>- beax with rate a (2.3) 

xaeb =>- xeab with rate f3 (2.4) 

baeb =>- beab with rate aj3 (2.5) 

ybey =>- yeby with rate 1 (2.6) 

abey =>- aeby with rate a (2.7) 
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ybea ==>- yeba wth rate j3 (2-8) 
abea ==>- aeba with rate a(3 (2-9) 

No other steps are allowed. In these equations the occupation of the four sites i, i+2 
is given and the dynamical step takes place between sites i and i + 1. We denote by x a 
site which is occupied by either an a-particle or by a hole e, and y represents a site which 
is either occupied by a 6-particle or by a hole e. The interaction between the two kinds of 
particles is attractive for a < 1 , (3 > 1 and is repulsive for a > 1 , (3 < 1. The other 
regions of the a, /3-plane describe dynamics with competing interactions. 

One can look at d > 1 dimensional systems as described by coupled chains which are 
directed along the driving field. The coupling between the chains allows for hopping of 
particles from one chain to another. As a result overtaking of a- and 5-particles in a single 
chain can take place by hopping via neighbouring chains. One may therefore model certain 
aspects of higher dimensional systems by studying a one-dimensional model with overtaking. 
The rate 7 of overtaking events is expected to be low. We thus consider a second version 
of the model where the following steps, involving overtaking, are permissible in addition to 
those given by Eqs. 



ea(e)be =^> eb(e)ae with rate 7 (2-10) 

aa(e)be ab(e)ae with rate /?7 (2-11) 

ba(e)be =^> bb(e)ae with rate 07 (2-12) 

ea(e)ba eb(e)aa with rate aj (2-13) 

aa(e)ba ab(e)aa with rate a/?7 (2-14) 

ba(e)ba bb(e)aa with rate ct 2 7 (2-15) 

ea(e)bb eb(e)ab with rate fij (2-16) 

aa(e)bb ab(e)ab with rate [3 2 ^ (2-17) 

ba(e)bb =^> bb(e)ab with rate a(3^ . (2-18) 

Another set of allowed processes are obtained from (|2.10|) -( 2.18D by interchanging a and b. 



The symbol (e) means, there can be a hole between the a- and the b- particle, which stays 
at the same position during the process. Note that the processes with a hole between the 
particles involve a next-nearest-neighbour interchange. 

III. MEAN FIELD APPROXIMATION 

In this section we discuss the stationary current in mean field approximation (MFA). 
The mean field approximation gives the same stationary current for both versions of the 
model. The stationary current of a particles in MFA can be written 

ja = Pa(l - Pbfil ~ P) + (a + P)PaPb(l ~ Pb)(l - P) + a[3p a pl{l - p) (3.1) 

The expression for the 6-particle current is similar. In terms of the ratio r introduced in Eq. 
(|2.1|) we get for the total current j = j a + jb 

j = p(l - p) + r(l - r)p 2 (l - p){2(a + [3 - 2) + (ot/3 - a - (5 + l)p} (3.2) 

The current exhibits a minimum in the conductivity if the relation 

2(a + f3-2) + (af3-a-f3+l)p<0 (3.3) 

is satisfied. 

Typical currents j as a function of r for different values of p are given in Fig.l. Although 
the current exhibits a minimum at r = 1/2 (which is most pronounced for a = (3 = 0), its 
value at the minimum is rather close to the current of the pure system ( j(l) = j(0) = p(l—p) 
), unlike the experimental results. In the next section we calculate the exact current for f3 — 
and demonstrate that the current at r = 1/2 can become much lower than its mean field 
value. 

IV. EXACT RESULTS AND SIMULATIONS 
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A. The model without overtaking 

We now discuss the first version of the model (rules ( |2.2|) - (|2.9|) ). Let us view the ring 
as a chain with sites 1, L and denote a configuration by n = {rii} where = a , b or e if 
site i is occupied by either an a-particle, 6-particle or a hole, respectively. Furthermore we 
define a frame F as the sequence of a- and 6-symbols which is obtained from a configuration 
n by removing all e-symbols. For example the frame associated with the configuration 
abebbeeaab... is abbbaab... . 

The dynamics of the first version of the model allow for jumps of a- or 6-particles into 
holes but do not allow for exchange of position of two particles. It may be viewed as the 
dynamics of holes moving to the left on a lattice certain sites of which are occupied by parti- 
cles. Consequently, the frame F_ changes in time only by means of cyclic permutations. The 
latter ones are possible because of the periodic boundary conditions assumed throughout 
this work. The dynamics are thus nonergodic, and the final state depends on the initial 
frame defined by the initial configuration rfi\ 

The case a = (3 = 1 . Although the experimentally observed anomaly is not expected 
here, a short discussion of this case seems instructive. There is neither repulsion nor at- 
traction between the a- and 6-particles for this choice of the parameters a and f3. The 
stationary probability distribution p s t(n) for this process is similar to the one for exclusive 
diffusion of one kind of particle which is well known [|20| . It assigns the same probability to 
all configurations with the initial frame F® or its cyclic permutations. The stationarity of 
this distribution can be easily checked by counting the number of incoming and the number 
of outgoing states for a given configuration. (An incoming state is a configuration which 
can change to the configuration of interest during one hopping event. Outgoing states are 
created by a single event taking place in the given configuration.) Note that the form of the 
stationary distribution depends strongly on the boundary conditions. The above discussion 
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is correct only for periodic boundary conditions. 

Obviously, since there is no interaction between the a- and 6-particles beyond the hard 
core term, the current j does not depend on the ratio r. As in the case of one kind of 
particles it is exactly the same as calculated in MFA: 

j(r,p)=p(l-p) . (4.1) 



The case (3 = (a ^ 0). This is a nontrivial but exactly soluble case, in which perfect 
repulsion between the a- and ft-particles takes place. Since f3 = 0, no a6-bonds (including 
ki-bonds) are created. But since a ^ afr-bonds may be broken. Hence the system runs 
into configurations for which the number of afe-bonds is minimal. The stationary state is 
thus characterized by the number p bo L of still existing afr-bonds. 

To study the stationary distribution in detail let us first define the quantity p^ as 
Pal = n ab/L where is the number of afr-bonds in the initial frame . (Note that 
due to the periodic boundary conditions we have to include in n^l bonds which may exist 
between site L and site 1.) The number of afr-bonds in the frame does not change in time 
for the only allowed changes of the frame are cyclic permutations. A configuration, however, 
has generally less afr-bonds than the frame, since holes may be located between the particles 
of such a bond. In order to allow for a configuration without any afr-bond the system needs 
at least as many holes as number of afr-bonds in the frame. Hence, if 1 — p < p®, the system 
runs into a configuration with a density p bo = p^ — 1 + p of a6-bonds and all holes are 
stuck between a- and 6-particles. Since any change of configuration takes place by means 
of a hopping hole, the configuration the system runs into does not change in time and has 
a vanishing current. On the other hand, if 1 — p > p^ b \ the number of holes exceeds the 
number of a6-bonds in the frame. The system evolves into configurations with no a6-bonds, 
i.e., where a- and 6-particles are separated by at least one hole. The number of holes free to 
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hop is given by rif = pjL with 

PS = 1 " P ~ P2 ( 4 - 2 ) 

It turns out that for 1 — p > p^ the distribution assigning the same probability to all config- 
urations with no afr-bonds and with the frame or its cyclic permutations is stationary. 
This can be seen by noting that the number of incoming states for any such configuration n 
is equal to the number of states to which the configuration can evolve. 

These simple considerations enable one to derive a general expression for the current 
j(r,p). Let us consider the case of nonvanishing current (1 — p > p^). A jump between two 
neighbouring sites, say site 1 and site 2, occurs with rate 1 if the left site is occupied by a 
particle and the right one by a free hole (a hole which is not stuck between an a- and an 
6-particle). Therefore j(r,p) equals the probability of finding ae/ or bef at sites 1 2, where 
e/ denotes a free hole. It reads 

j(r, p) = prob (ni = a or n\ = b) prob (n 2 = ef\ri\ = a or n\ = b) . (4.3) 

Here, the second term in the right hand side of this equation denotes the conditional 
probability of finding a free hole at site 2 given that site 1 is occupied by a particle. 
Due to translational invariance one has prob {n\ = a or ri\ = b) = p. The probabil- 
ity prob (n 2 = e.j\n\ = a or n\ = b) is given by the ratio of the number of configu- 
rations which have a free hole at site 2 and a particle at site 1 to the total number of 
configurations with a particle at site 1. The configurations contributing to these num- 
bers may have different frames. However, all these frames can be obtained from by 
cyclic permutation. Thus they all have the same number of afe-bonds, and we may re- 
strict ourselves to configurations with the frame The above probability is given by 
prob {ri2 — ef\ni — a or n\ — b) — X/Y, where Y is the number of configurations with the 
frame and a particle at site 1, and X is the number of those configurations which in 
addition have a free hole at site 2. All configurations to be counted may be constructed by 
first inserting a hole in between any pair of a6-particles in the frame . One then has to 
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distribute the remaining n/ = p/L free holes in between the particles in a way that site 1 is oc- 
cupied by a particle. Thus Y is equal to the number of ways of distributing nf indistinguish- 

i n+nf — l i 

able holes in n = n a + states. It is given by Y = \ \ . Out of these configurations 

n+rif— 2 

X = Y — have a free hole at site 2. The second term on the right hand side of this 



f 

equation gives the number of configurations with no free hole at site 2. The current j(r,p) 

for 1 — p > is therefore given by: j(r, p) — p { f f I — I / I } I 1 

V n f J \ n f J \ n f 
Simplifying the binomials results in the following expression for the current: 



. + PPf r ! for 1 - p > p { 2 
j(r,p) = { (4.4) 



'ah 

otherwise 

,(*) 



The density of free holes, pf, is a function of p and of p a l , as expressed in (|4.2j ). In the 
thermodynamic limit the L _1 -term in the expression for the current vanishes. 

Consider now the case of random initial conditions in which the initial configuration 
is created by uniformly distributing p a L a-particles and pt,L 6-particles on a lattice of length 
L. The average number of afr-bonds in the frame F} 1 ' corresponding to this initial condition 
is 2r(l — r)pL. For L —>■ oo the relative fluctuations of this quantity vanish. Consequently, 
we find for any initial configuration with a- and b-particles uniformly scattered on a large 
lattice: 

pg = 2r(l-r)p . (4.5) 
Combining this with Eq. (|4.2|) and (f4.4|) we obtain: 



P l ^ilt{ )] p < Pc (r) 



j(r,p) = {' L -^ l - r) (4.6) 
otherwise 

where 

p c {r) = 1/[1 + 2r(l - r)] (4.7) 

is the critical density. The expressions for the current j and the critical density p c are 
independent of a (as long as a ^ 0). Fig.l shows the current as a function of r for different 
densities p. For densities p < p c (r = 1/2) = 2/3 the current is nonzero for any value r. 
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This result is rather different from the mean field current ( |3.2| ). While the latter is a 
function of a the exact expression for the current j is independent of a. Fig. 1 compares the 
exact result with the mean field current ( |3.2| ) for a = 0. This value is chosen since it corre- 
sponds to the most pronounced minimum of the conductivity in mean field approximation. 
The discrepancy between mean field and exact results indicates that the deep minimum is 
related to correlations. 

The case (3 > (a ^ 0) . For nonvanishing (3 the stationary distribution is more difficult 
to calculate and one has to resort to numerical simulations of the model. However, a few 
properties of the j(r, p)-function can be seen easily. For r = and r = 1 the current is 
equal to the mean field current given in Eq. (|4.1| ) since in both cases only one kind of 
particles performs exclusive diffusion. The stationary solution for that process is known to 
be uncorrelated ||. Furthermore the j(r)-function is symmetric with respect to reflection 
about the r = 1/2-line because the model is defined in way that a- and b-particles play the 
same role. 

The stationary current for (3 > does not vanish for high densities, as in the case (3 = 
since there are always possible hopping events. Consequently, we expect a qualitatively 
different behaviour for p > p c as compared to the (3 = 0-case. Replacing (3 = by a finite (3 
the current is increased for any value of p and r. 

Computer simulations of the process were performed by letting particles hop stochasti- 
cally on a lattice of L = 1000 sites. Averages are calculated as time averages for a particular 
realization as is done in experiments. Besides, time averages and ensemble averages coincide 
if the stationary probability distribution is chosen properly, i.e. in the subspace of the phase 
space which is actually reached by the system. 

We have carried out simulations for the case < (3 « 1 and p = 0.7 > p c which seems 
to reflect the experimental results for the current ||] quite well. The parameter a is chosen 
to be 1 but it is supposed to play a minor role as long as a ^ 0. Results are shown in 
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Fig. 2. Tuning (3, which depends on the temperature via (3 = exp(—5E/kBT), any ratio 
j(r = 0)/j(r = 0.5) can be obtained for p > p c . Here 8E is the energy barrier one has to 
overcome by creating an afe-bond. The model is therefore capable of exhibiting the deep 
minimum of the conductivity which is observed experimentally. 



B. The model with overtaking 



We next study the second version of the model (rules @)-(|jJ)) where 7^0. bmce 



overtaking is allowed in these dynamics, the frame undergoes noncyclic permutations of 
particles unlike in the 7 = case. 

Let us discuss the case of perfect repulsion (/3 = 0) in detail. By simple state-counting 
it can be shown that a distribution assigning the same probability to any occurring con- 
figuration is stationary. But while only configurations with the initial frame Fy' and its 
cyclic permutations are allowed for the first version of the model, here there is a broader 
distribution of frames which are reached by the dynamics. As in the previous version, for 
sufficiently large density p the model is expected to run into a configuration in which all 
holes are stuck between a- and 6-particles. This state has a vanishing current. On the other 
hand, for densities lower than some critical density p c (r) the number of holes exceeds the 
number of afr-bonds in the initial frame. The system therefore has free holes which generate 
nonvanishing current. 

For densities satisfying p < p c {r) rules (|2.10|) - (|2.18|) allow for changes in the frame via 



interchanges of a- and 6-particles. Any frame with a number of bonds less than the number 
of holes can be created during the dynamics. The stationary probability distribution p st 
assigns the same probability to any configuration without a6-bonds. 

To calculate the current associated with the stationary distribution we first consider the 
current averaged over all configurations whose frames have n a b = p a bL afr-bonds. As it was 
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shown in Section IV A the current averaged over all configurations with p and p ab is given 
by 

JWH^T^T < 4 - 8 > 

where 

Pf = 1 " P ~ Pab (4-9) 

is the density of the free holes. The steady state current is thus obtained by averaging 
j(p, Pab) over all possible p a b or n ab . Since all allowed configurations have the same weight one 
has to find the number of configurations corresponding to n ab . To this end we consider first 
the number of frames associated n ab . We then calculate the number of ways of distributing 
rif = pfL free holes on a given frame, and obtain the probability of having a frame with 
n a b afr-bonds. Note that, as explained at the beginning of Section IV A, a frame is defined 
by taking a configuration n = {rii, ...,n^} and removing all holes. This defines a sequence 
of A intervals of a-particles alternating with B intervals of 6-particles. The first and the 
last interval may either be of the same or of different type. If the intervals at both ends 
of the lattice are the same, say a, type then A = n ab /2 + 1 and B = n ab j2. Similarly if 
both end intervals are of 6-type then A = n ab /2 and B = n ab /2 + 1. On the other hand if 
the two end intervals are of different types then A — B — n ab /2. The number of ways of 

(rta-l \ 
, and similarly the number 
A-lJ 

of possibilities of arranging the n b = (1 — r)pL 5-particles in B groups is I & ] . Therefore 

V"-'/ 

in each of the above cases the number of possible frames is I ] ( b ] . We now have 

\a-iJ \b-iJ 

to find the number of ways of distributing rif = pfL free holes in each frame. In the case 
where both ends of the frame are of the same type this number is given by while 



n 



f 



, . i L—riab \ I L—n a if—l . 

in the case where the two ends are different it is given by [ + ]■ Thus 



the statistical weight associated with n ab takes the form 
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ab) 



n a -l J / n 6 -l + n a -l \ / n b -l \ I L-n ab 
nab/'ij VW2-1/ \n ai ,/2-l/ \n ai) /2/ \ 

+2 I ™ a 1 | ( nfc_1 | [ | L ~ riab | _|_ / i_ "ab-i 

yn at ,/2-iy \n ab /2-l) \ n f J \ n f 

2L I L-n ab -l \ I r P L-X \ ( (l-r)pL-l 
n ab \{l-p)L- nab ) \n ab /2-l) I n ab /2-l 



(4.10) 



Using Eqs. (]4.8|)-(3~TD|) we find for the stationary current: 



pz- 1 ET=oT^mjL^f( 2k ) fOT P<Pc(r) 







otherwise 



(4.11) 



where Z = Y^kLo f {2k) is a normalization constant. 

Equations ( f4.10|) - ([4.11|) give the stationary current for a lattice of arbitrary length L. In 
the thermodynamic limit the weight function f{n a b) is sharply peaked around n a b = p a b*L 
where p a b* is a solution of the following equation: 



= r(l - r)p(l - pf - -(1 - p)[l - p + 4r(l - r)p] p ab * 



J r[^-\p J rr{l-r)p] {p. 



*\2 



This yields p a b* is a function of p and r(l — r) which has its maximum at r 
given value of p. For example , solving Eq.( 4.12 ) for r = 1/2 we find 



(4.12) 



1/2 for any 



1 



pab (r 



P) 



1 1 



l-2p(l-p) . 



(4.13) 



Hence, for an infinite lattice the number of bonds in the frame is 
current is 



i,*L and the stationary 



j{r, p) 



P 




X-p-p ab *{r,p) 
l-Pab*{r,p) 



for p < p c {r) 
otherwise 



(4.14) 



Let us now consider the random initial conditions defined in Section IV A and try 
to estimate the critical density p c {r) below which the current of the system is nonzero. 
Clearly, if the concentration p^ = 2r(l — r)p of a6-bonds in the initial frame satisfies 
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p^ < 1 — p., the system has free holes which generate ergodic dynamics leading to a current 



j as given by Eqs. ( 4.13|) and (f4.14 ). This yields a first lower bound for the critical density 



p c (r) > 1/[1 + 2r(l — r)]. However, due to the particle exchange mechanism which exists 
for 7 > the system may exhibit a finite current even for p^ ^ 1 — p. This may easily 
be seen by considering the limit 7 C 1. In this limit and for time scales shorter than I/7, 
basically all initially free holes are caught by a6-bonds, and thus stop moving. However for 
longer time scales, where particle exchange processes take place, some of the holes which 
are stuck in afe-bonds are released. For example a sequence abab may evolve into aabb thus 
reducing the number of a6-bonds by two and creating free holes. And a sequence aabbaabb 
may in principle evolve into aaaabbbb, but this requires the existence of free holes to begin 
with. If the only holes in this sequence are those stuck in the afe-bonds, particle exchange 
processes do not take place. However, if in addition there are some free holes, the sequence 
may change. If we take into account only changes which do not require the existence of free 
holes the number of a6-bonds in the frame is reduced and becomes p ab L with 

At = Pal ~ [4r 2 (l - rf + o(r 3 (l - rf)]p . (4.15) 

We thus expect that for p ab < 1 — p the system exhibits a nonvanishing current. This yields 
the following lower bound for p c (r): 

p c (r) > 1/[1 + 2r(l -r) -4r 2 (l - r) 2 ] + o(r 3 (l - rf) . (4.16) 

The current j(r,p) is shown in Fig.3 as a function of r for various values of p. In the 
figure p c {r) is determined by Eq. ( |4.16| ) to second order in r(l — r). 

The dynamics of the hopping holes have a time scale of order 1 whereas the reordering of 
the frame has a characteristic time much larger than one since we consider the case 7 -C 1. 
Hence we observe the following scenario: The system runs into the stationary state of the 
first version of the model (see Fig. 1) which decays very slowly to the final state (Fig. 3). 
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V. CONCLUSION 



A simple model describing the transport properties of mixed ionic conductors has been 
introduced and analyzed. The model exhibits a minimum conductivity for equal concentra- 
tions of the two species of particles (r = 1/2), a result which is compatible with experimental 
observations in mixed alkali glasses. A mean field approximation yields the correct qual- 
itative behaviour of the conductivity, but it fails to explain the very low conductivity for 
r = 1/2. It has been demonstrated by exact solution and numerical simulations in the 
strong repulsion limit that the conductivity corresponding to the model is indeed very small 
at r = 1/2 in accordance with experimental observations. 

The model can be extended to study the case in which particles move into both directions 
with arbitrary rates. The exact results obtained above for the case (3 = are easy to gener- 
alize: The probability distribution which assigns the same probability to all configurations 
with no bonds is stationary even when hops in both direction take place. Moreover, the 
rates p and q of hopping to the right and left, respectively, may be taken as time dependent 
without changing the stationary state. The current is then \p(t) — q{t)]j, where j is given 
by Eq.( [4.4|) or Eq. (|4.6| ) for 7 = and by Eq. (|4.11| ) or Eq. (|4.14|) when overtaking is included, 
i.e. for 7 7^ 0. Obviously, the conductivity does not depend on the frequency if the driving 
field is harmonic. 



Acknowledgments 

One of us (S. Sandow) gratefully acknowledges financial support by the Minerva Foundation. 



14 



[1] A comprehensive review of early literature, see D.E. Day J. Non-Cryst. Solids 21, 1976 

[2] P. Maas, A. Bunde and M.D. Ingram Phys. Rev. Lett. 20 3064, 1992 

[3] A. Bunde, M.D. Ingram, P. Maass and K.L. Ngai J. Phys. A: Math. Gen. 24 L 881, 1991 

[4] A. Bunde, P. Maass and M.D. Ingram Ber. Bunsenges. Phys. Chem. 95 977, 1991 

[5] H. Spohn 1991, Large Scale Dynamics of Interacting Particles, 1991 

[6] S. Katz, J.L. Lebowitz and H. Spohn, J. Stat. Phys. 34 497, 1984 

[7] P.M. Richards, Phys. Rev. B 16 1393, 1977 

[8] T.M. Liggett Interacting Particle Systems (Berlin: Springer), 1985 

[9] P. Meakin, P. Ramanlal, L.M. Sander and R.C. Ball, Phys. Rev. A 74 5091, 1986 

[10] J. Krug and H. Spohn in Solids Far From Equilibrium: Kinetic Roughening of Growing Surfaces 
ed C Godreche (Cambridge: Cambridge University Press), 1991 

[11] D. Kandel and D. Mukamel, Europhys. Lett. 20 325, 1992 

[12] L.-H. Gwa and H. Spohn, Phys. Rev. A 46 844, 1992 

[13] B. Derrida, E. Domany and D. Mukamel, J. Stat. Phys. 69 667, 1992 

[14] B. Derrida and M.R. Evans, J. Physique 13 311,1993 

[15] G. Schiitz and E. Domany, J. Stat. Phys. 72 277,1993 

[16] S. Sandow and G. Schiitz, Europhys. Lett., 26 7, 1994 

[17] B. Derrida, M.R. Evans, V. Hakim and V. Pasquier, J. Phys. A: Math. Gen. 26 1493, 1993 

[18] G. Schiitz, J. Stat. Phys. 71 471,1993 

[19] B. Derrida, M.R. Evans and D. Mukamel, J. Phys. A: Math. Gen. 26 4911, 1993 

[20] F. Spitzer, Adv. Math. 5 246, 1970 

15 



[21] B. Derrida, S.A. Janowsky, J.L. Lebowitz and E.R. Speer E R, Europhys. Lett. 22 651,1993 



Captions to the figures 

Fig.l: Stationary current j for the model without overtaking as a function of the ratio 
r for (3 = and different values of the density p, solid lines: exact results (a arbitrary), 
dashed lines: mean field approximation (a = 0); (1): p = 0.15; (2): p = 2/3 ; (3): p = 0.8 

Fig. 2: Stationary current j for the model without overtaking as a function of the ratio r for 
p — 0.7 , a — 1 and different (3 l s obtained by simulations (The labels are the values of (3, 
and the lines are interpolations.) 

Fig. 3: Stationary current j for the model with overtaking as a function of the ratio r for 
(5 — 0, arbitrary a and different values of the density p, exact results 
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